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We extend the theory of Kimberly and Magueijo for the spacetime variation of the electroweak 
couphngs in the unified Glashow-Salam- Weinberg model of the electroweak interaction to include 
quantum corrections. We derive the efi'ective quantum-corrected dilaton evolution equations in the 
presence of a background cosmological matter density that is composed of weakly interacting and 
non-weakly-interacting non-relativistic dark-matter components. 

I. INTRODUCTION 



The studies of relativistic fine structure in the absorption lines of dust clouds around quasars by Webb et al., [1, 2, 3], 
have led to widespread theoretical interest in the question of whether the fine structure constant, aem = /he, has 
varied in time and, if so, how to accommodate such a variation by a minimal perturbation of existing theories of 
electromagnetism. These astronomical studies have proved to be more sensitive than laboratory probes of the constancy 
of the fine structure 'constant', which currently give bounds on the time variation of aem/ctem = — 0.4±16x lO^^^yr"^, 
[4], \aemloiem\ < 1-2 X 10~^^?;r~^, [5], diem/oiem = — 0.9 ± 2.9 X 10~^^yr~^ , [6] by comparing atomic clock standards 
based on different sensitive hyperfine transition frequencies, and ctem/ctem = —0.3 ±2.0 x 10^^^yr~^ from comparing 
two standards derived from 1S-2S transitions in atomic hydrogen after an interval of 2.8 years [7]. The quasar 
data analysed in refs. [1, 2, 3] consists of three separate samples of Keck-Hires observations which combine to 
give a data set of 128 objects at redshifts 0.5 < z < 3. The many- mult iplet technique finds that their absorption 
spectra are consistent with a shift in the value of the fine structure constant between these redshifts and the present of 
Aaem/c^em = [cemiz) — aem] / o^em = — 0.57±0.10x 10~^, whcrc ttem = c«em(0) is the prcscut value of the fine structure 
constant [1, 2, 3]. Extensive analysis has yet to find a selection effect that can explain the sense and magnitude of the 
relativistic line-shifts underpinning these deductions. Further observational studies have been published in refs. [8, 9] 
using a different but smaller data set of 23 absorption systems in front of 23 VLT-UVES quasars at 0.4 < z < 2.3 
and have been analysed using an approximate form of the many- mult iplet analysis techniques introduced in refs. 
[1, 2, 3]. They obtained Aaem / ctem = —0.6 ± 0.6 x 10"^; a figure that disagrees with the results of refs. [1, 2, 3]. 
However, reanalysis is needed in order to understand the accuracy being claimed. Other observational studies of lower 
sensitivity have also been made using OIII emission lines of galaxies and quasars. The analysis of data sets of 42 and 
165 quasars from the SDSS gave the constraints Aaem/ctem = 0.51 ± 1.26 x 10~^ and Aaem/oiem = 1-2 ± 0.7 x lO"'^ 
respectively for objects in the redshift range 0.16 < z < 0.8 [10]. Observations of a single quasar absorption system 
at z = 1.15 by Quast et al [11] gave Aaem/ckem = —0.1 ± 1.7 x 10"® , and observations of an absorption system at 
z — 1.839 by Levashov et al [12] gave Aaem/ceem = 4.3 ± 7.8 x 10"®. A preliminary analysis of constraints derived 
from the study of the OH microwave transition from a quasar at z = 0.2467, a method proposed by Darling [13], has 
given Aaem/aem = 0.51 ± 1.26 x 10"'', [14].A comparison of redshifts measured using molecules and atomic hydrogen 
in two cloud systems by Drinkwater et al [15] at z = 0.25 and z = 0.68 gave a bound of Aaem/ctem < 5 x 10~® 
and an upper bound on spatial variations of daem/c(em < 3 x 10"® over 3 Gpc at these redshifts. A new study 
comparing UV absorption redshifted into the optical with redshifted 21cm absorption lines from the same cloud in 
a sample of 8 quasars by Tzanavaris et al [16]. This comparison probes the constancy of a^gpTne/rrip and gives 
Aaem/ctem = 0.18 ± 0.55 X 10"^ if we assume that the electron-proton mass ratio and proton g-factor, gp, are both 
constant. 

Observational bounds derived from the microwave background radiation structure [17] and Big Bang nucleosynthesis 
[18] are not competitive at present (giving Aaem/ctem ^ 10~^ at best at z '--^ lO'^ and z ~ 10® — 10^°) with those 
derived from quasar studies, although they probe much higher redshifts. 

Other bounds on the possible variation of the fine structure constant have been derived from geochemical studies, 
although they are subject to awkward environmental uncertainties. The resonant capture cross-section for thermal 
neutrons by samarium-149 about two billion years ago (z ~ 0.15) in the Oklo natural nuclear reactor has created a 
samarium-149:samarium-147 ratio at the reactor site that is depleted by the capture process ^^^Sm + n — s-^^" Srn + ^ 
to an observed value of only about 0.02 compared to the value of about 0.9 found in normal samples of samarium. 
The need for this capture resonance to be in place two billion years ago at an energy level within about QOmeV of its 
current value leads to very strong bounds on all interaction coupling constants that contribute to the energy level, as 
first noticed by Shlyakhter [19, 20]. The latest analyses by Fujii et al [21] allow two solutions (one consistent with no 
variation the other with a variation) because of the double- valued form of the capture cross-section's response to small 
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changes in the resonance energy over the range of possible reactor temperatures: /^a^m/ c^em = —0.8 ± 1.0 x 10~^ or 
Aofem/o^em = 8.8 ± 0.7 X 10~^. The latter possibility does not include zero but might be excluded by further studies 
of other reactor abundances. Subsequently, Lamoureax [22] has argued that a better (non-Maxwellian) assumption 
about the thermal neutron spectrum in the reactor leads to 6(7 lower bound on the variation of Aaem/ctem >4.5xl0~* 
at 2 ~ 0.15. 

Studies of the effects of varying a fine structure constant on the /3-decay lifetime was first considered by Peebles 
and Dicke [23] as a means of constraining allowed variations in by studying the ratio of rhenium to osmium in 
meteorites. The /3-decay \^ Re Os + + e~ is very sensitive to agm and the analysis of new meteoritic data 

together with new laboratory measurements of the decay rates of long-lived beta isotopes has led to a time-averaged 
limit of Aaem/ctem = 8 ± 16 X 10"'' [24] for a sample that spans the age of the solar system {z < 0.45). Both the 
Oklo and meteoritic bounds are complicated by the possibility of simultaneous variations of other constants which 
contribute to the energy levels and decay rates; for reviews see rcfs. [25, 26]. They also apply to environments within 
virialised structures that do not take part in the Hubble expansion of the universe and so it is not advisable to use 
them in conjunction with astronomical information from quasars without a theory that links the values of aem in the 
two different environments that differ in density by a factor of O(10^°). 

In order to interpret these investigations it is essential to be in possession of a self-consistent theory of the spacetime 
variation of aem, analogous to the Brans-Dicke theory [27] for the variation of G, from which to derive further 
observational consequences of any inferred variation. In the past, in the absence of any such theory, there has been 
a tendency to produce limits on variations of couplings other than agm by simply 'writing in' a time variation of the 
coupling into the equations that hold when it is constant. Many of the quoted experimental bounds on the variation 
of non-gravitational 'constants' that appear in the literature have been arrived at in this way. Also questionable is the 
use of laboratory bounds to limit possible variations of 'constants' on extragalactic scales without any theory of the 
link between the two domains of variation. Detailed discussions of this problem when G and a vary have been made 
in refs. [28, 29]. A further problem is the likelihood that if one coupling constant, like aem, varies then others will 
vary also, especially in the presence of any pattern of unification [15, 30]. Most deductions of bounds on constants 
assume that only one constant is varying. 

In order to deal with these problems of simultaneous variation and spatial variation consistently, it is necessary 
to have self-consistent theories in which 'constants' vary through the dynamics of scalar fields, which gravitate and 
conserve energy and momentum. Sandvik, Barrow and Magueijo (SBM)[31] have extended Bekenstein's (B) [32] 
original generalisation of Maxwell's theory to include general relativity . This resulting BSBM theory provides a 
framework for studying varying aem, said can be easily generalised to study the simultaneous variation of aem and 
G, as in ref. [33]. A number of detailed cosmological studies of the behaviour of this theory have been made in 
refs. [31, 34, 35, 36, 37, 38]. A further step has recently been taken by Kimberly and Magueijo [39] who have 
extended the BSBM theory to create a generalisation of the Glashow-Salam- Weinberg theory to allow variation of the 
electromagnetic and weak couplings. This allows the consequences of electroweak unification to be investigated self 
consistently for the first time. This approach has also been applied to the strong interaction alone by Chamoun et al 
[40] and it is likely that a further step could be taken to investigate the consequences of the simultaneous variation of 
all gauge couplings in a grand unified theory. In this paper we determine the quantum corrections to the Kimberly- 
Magucijo unfied electroweak models and formulate the propagation equations for the two scalar fields that carry the 
allowed time variations in the electromagnetic and weak couplings in a universe populated by the particle spectrum 
of the standard model. 



II. THE SIMPLEST MODEL 

There has been strong interest in a class of scalar theories where Uem = /he can vary [31, 32, 34, 35, 36, 37, 38, 
41, 42, 43, 44]. These theories reduce to Maxwell's equations and general relativity in the limiting case of no variation 
in the fine structure constant. In order to evaluate the constraints introduced by a programme of unification, it is 
important to extend these simple models to include electroweak and grand unification. Kimberly and Magueijo [39] 
have proposed an extensions to the Glashow-Salam- Weinberg (GSW) electroweak theory in which the weak couplings 
can also vary in space and time and which reduces to the standard GSW theory in the limit of no variation. 

A. A Single-Dilaton Theory 

The first model (KM-I) contains a single dilaton and allows both aw '■= Qw ^'^d aem to vary but their ratio, 
and hence the mixing angle 6w, are true constants. In this model the electroweak sector of the Standard Model is 
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described by the following lagrangian density: 
where 



Vf^„ = 2wiy^^^- g^^N^,^■w^, (2) 
y^i- = 2y[^,^], (3) 

[di^ - ^9w^ ■ - ^^gYy^^^ * (4) 

and $ is the Higgs field. The dilaton field is (p and w is a dimensional parameter with units of (mass)'^. Since this 

theory is pcrturbatively non-renormalisable, wc would like lo ^ O (^Afp/^ so that the dilaton enters at the same level 

as gravity, and we are justified in ignoring any quantum fluctuations of the dilaton field when quantising with respect 
to the other fields. The auxiliary gauge fields, and y^, take values in the adjoint representations of 5u(2) and u(l) 
respectively. They are not the physical gauge fields, which will be denoted by capital letters, but are related to them 
by the transformations 

gw^n = gw'Wfi, (5) 
gvyi^ = qyY^, (6) 
gw = gwe"^, (7) 
QY = gYe'P. (8) 

The distinction between the physical and auxiliary fields will only become important at the quantum level (see 
section 2 below). When written in terms of these auxiliary gauge fields the covariant derivatives which act upon 
matter species are independent of This makes it simpler to derive the classical field equations. The physical gauge 
couplings, gw and qy, are dynamical whereas the auxiliary couplings, gw and gY, are arbitrary constants. At tree 
level the Fermi constant, Gp, and the fermion masses do not vary, whereas the W and Z boson masses do. 

We will also define the physical field strength tensors by 

W^^ := w^^, (9) 

gw 

Y^. := — y^.. (10) 

gY 

These field strengths reduce to the standard definitions of the weak and hypercharge field strengths with gauge 
couplings gw and gY respectively whenever the dilaton field ^ is constant. 



B. A Two-Dilaton Theory 



In the second model (KM-II)proposed in ref. [39] a second dilaton field is added. This results is the weak mixing 
angle, Qw i becoming a dynamical quantity. The lagrangian density of the elcctroweak sector in this model is: 



KM-II 



ljJ2 



(11) 



The definitions (2)-(4) still hold, as do relationships (5)-(6) but the physical coupling constants are now related to 
their auxiliary values by the transformations 



gw = gwe 
gY = SYe^, 

tan^vv' 



- — := ta,now 

gw 



(12) 
(13) 

(14) 



The dilaton fields are ip and x and their respective dimensionful scales are wi and W2- As remarked above, we would 
like - O (M|,;). In 
masses are dynamical. 



like cjj ^ O (^MpiY In accord with KM-I, Gp and the fermion masses are constant at tree-level whereas the boson 
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C. Symmetry Breaking 

At low energies and temperatures, the most important feature of the GSW electroweak model is that the SU{2) l x 
U{\)y symmetry of the lagrangian is spontaneously broken to ?7(l)em via the Higgs doublet, ^, assuming a vacuum 
expectation value, which minimises its potential. At tree- level this value is 

*o = . (15) 
A perturbative expansion about this vacuum can be written as 

*=(. + %i)^ (16) 

Expanding out the kinetic Higgs term gives: 

(D^*) = \h,H'^ + ^-^^ [(Wl)2 + (W^)^] + !^Cg^^l-g^y^f. (17) 

When written in terms of the physical gauge fields, the broken-phase boson fields of both KM-I and KM-II are given 
by the usual formulae: 

W^:=^{Wl±iWl), (18) 
_ 9wWI-9yY, 
V9w+9y 



V9w + 9y 



A, := ^ (20) 



Hence, the tree-level boson masses and their dilaton field dependence can read off as(17): 

(21) 



Mw = 




Mz = 


-^\/9w + 9y 


Ma = 


0. 



(22) 
(23) 



D. Classical Field Equations 

In order to make and test the predictions of these theories we need to know the field equations. In both KM-I and 
KM-II, the Einstein and Yang-Mills equations are: 



rpmatter i rpdilaton 1 

1 9w 9y J 



(24) 



where and TJ^^, are the standard Yang- Mills energy-momentum tensors written in terms of the auxiliary fields and 
couplings. In KM-I the dilaton conservation equation is 

n^ = - (w^. • w**^ + y^^y^n = (W^. • + Y^^Y^^ . (27) 
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and in KM-II we have conservation equations for the two fields: 

-:^^e-2vw^.•w'''' = --^W^,•W^^ (28) 

The conservation equations for other matter fields, like perfect fluids, are unchanged. Although this system represents 

the full classical field equations, their current form is not very useful. In order to do cosmology, and make experi- 
mentally testable predictions, we need to understand how the right-hand sides of (27, 28, 29) depend on macroscopic 
quantities such as the background matter density and the value of the dilaton field. It is this question that we address 
in the next section. 

III. THE DILATON-TO-MATTER COUPLING 

Most tests which we might wish to apply to these, and other, varying-constant theories require knowledge of how 
the dilaton fields will evolve in the presence of some background matter density, p, and pressure, P: 

p ^^^O(matier) ^ p ^ i{matter) 

where ( • ) denotes the quantum expectation. In order to understand the dilaton evolution, we must therefore evaluate 
the terms on the right-hand sides of (27)-(29) under the( • ) operation. Those terms all consist of terms quadratic in 
the Yang-Mills field strengths, i.e. W^y • W"" and F^,,F'"', henceforth we refer to them collectively as F^^F'^'' or F^ 
terms. 

When we quantise this theory (leaving (/? as a classical field) we must do so with respect to the physical gauge fields 
rather than the auxiliary ones, since it is only the physical fields whose kinetic terms possess the correct normalisation. 
This feature is important if we want the renormalisation procedure to go through as usual when (p = const. It is for 
this reason that we have labelled the capitalised fields as physical. 

A. Derivation 

We need to understand how these (F"^) terms depend, in a functional sense, on the dilaton fields and the background 

matter density, p. For simplicity, we consider the case of a single Dirac fermion, ijj, of mass m, coupled to a U{1) 
gauge field with Higgs scalar $ and dilaton (p. The effective lagrangian, including 't Hooft's gauge-fixing term, is: 

£„,o = -IF^^'F^^ + II {d^ + i^A^) $f - A ($t$ _ (31) 
£/ix = -^(S''A^ + eMm($))' (32) 

C^fi = -0 {il^di, - to) V (33) 
Cint = -e (v) i^l^Mt. (34) 

The ghost fields have been excluded since in this gauge they decouple from the abelian sector. By taking a spacetime 
region TZ, of volume V and temporal extent T, that is large compared to the scale of quantum fluctuations of the 
matter and gauge flelds, and yet small when compared to the scale over which the dilaton field varies, we can quantise 
this theory inside TZ by taking ip = const and defining the partition function in the usual manner. For energies far 
below the Higgs mass, niH, we can ignore quantum Higgs fluctuations in the first approximation. Doing this and 
integrating out the gauge fields we determine that 

-i {F^.F^ =^ Iff d^xd^yd^z(dlV3l^s%x - y) (i^Yi^) {y)V%^-^<>- {x - z) (i^Yi^) (z) 

(35) 
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where the gauge-field propagator is given by: 



(1-0 fc'^fc^ 



fc2 - 

The quantum expectation operator, (•)^, used in equation (35), is defined thus: 

/ [d^] [dtp] A (Vj, iJj) cxp / / d^ard^yJ'' (a;) (a; - y) J'' {yj^e'I ^*'''^*-° 

(A (V', V') := J [d^] [d^] exp (-if / / d^xd^yj- (:c) (a; - y) J'' (y))e'/ d%/:^,o ^^^^ 

The / [dtp] represents a functional integral with respect to the f/'-field, and we have defined J'^{x) := tp{x)j'^tp{x). The 
term with an exponent that is quadratic in the J's in equation (36) encodes the self-energy corrections to the fermion 
and gauge boson propagators. By writing both the fermion and gauge boson propagators as full propagators we are 
justified in dropping the aforementioned exponential term because its effect is of sub-leading order. Finally then, by 
writing this expression in momentum space we reduce it to the more transparent and manifestly gauge invariant form: 



-i(i^,.Fn = ^/ ^p^[(p^-Ml(^)),''^+n7(p,e^(^)])-^(>(p)J''(-p))^^, (37) 

where is the momentum space representation of and U.'^ is the vacuum polarisation of the gauge boson. The 
expectation here is defined with respect to the partition function 

Z^^ := J [d^j] [dV]exp(^i J d''xd''yij{x)K{x-y)tP{y)^ 

with K{x — y) denoting the inverse of the full fermion propagator. 



B. Interpretation 



By Wick's theorem it is apparent that the remaining quantum expectation in (37) contains two distinct contribu- 
tions. The first contribution will always involve boson exchange between two fermionic particles and so is proportional 
to (p — 3P)^ /m^ to leading order. The second contribution results from a single fermionic particle admitting and 
reabsorbing a gauge boson (i.e. the process that results in the fermion's self-energy). This second term will clearly be 
proportional to m? {p — 3P) to leading order. In almost all cases of experimental interest {p — 3P) jm^ <^ 1. Only in 
objects whose density approaches that of nuclear matter does it fail to hold and this naive quantisation procedure is 
not suitable for dealing with such high-density backgrounds. For this reason we drop the contribution due to photon 
exchange. 

When the perturbation theory holds it is appropriate to expand (37) as a series in the gauge coupling, {(p) , giving 



-l{F^^Fn = e\i^)ip-SP){l + 0{e'lne)){l + 0^^)) (38) 
where, for this model, C is a defined by: 



d^p d*q ,/jy9-7 + 'Ti {p - q) • 1 + 



m 



C(^)(p-3P) = / / -^-^ o^Al 9 7. . J (39) 

When inter-fermion strong interactions are introduced, confinement may occur. If this happens the trace factor in 
(39) will take on a more complicated form. At leading order, however, it will remain independent of so long as 
^QCD is. At densities much lower than the nucleon mass, the right-hand side will be proportional to the hadronic 
energy density. 



C. The Unbroken Symmetry Case: Ma = 



It is helpful to have a physical interpretation of the C parameter. When dynamical symmetry breaking docs not 
occur (and so we have M\ = 0) the right-hand side of (39) is proportional to the 1-loop fermion field self-energy 
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resulting from its interaction with the A^^ gauge field. Defining this self-energy to be Sm {ip) we have: 

Sm^ (ip) Sfn^ , 

C = 2( V 2 = = (40) 

In this case C is ^p- independent at leading order. 5fh is defined as the electromagnetic mass correction when the electric 
charge is some value e = const. 



D. The Broken Symmetry Case: Ma # 

When M\ ^ the above interpretation of C, in terms of the self-energy mass-correction will, in general, fail. The 
correct physical interpretation will depend heavily on the size of jgj^- We consider the three cases: 

• <C 1. The dominant contribution to the Q integral will come from momenta ^ M\, and so we can ignore 
at leading order and interpret C just as we did in the Ma = case. Hence, C is dilaton independent to 

leading order. 

• ^ 1. The dominant contribution to the (39) will come from momenta m'^. Hence 

where Co is of the order the value of ( which we would have had if Ma = 0. In this case ( will vary with like 
1 

• = O (1). This is by far the most difficult case to analyse. In general ,C will have a leading-order dependence 
on the dilaton field through Ma (<^), but the precise form ( {(p) will depend on the nature of the trace term in 
equation (39). This in turn rests on the precise details of the microscopic physical model for the matter fields 
in question. This chain of complications means that a general prescription for the dilaton dependence of C, in 
this case, is not possible. 



IV. APPLICATIONS TO THE KIMBERLY-MAGUEIJO MODELS 



A. Ultra-Relativistic Matter 



At the very high energies required to restore the broken gauge symmetry in these theories, most matter species 
will be ultra-relativistic and we will assume this to be the case for all species. We will also assume that there is a 
discrete spectrum of spin states. Matter will therefore behave like black-body radiation with p ~ 3P. It is clear from 
equation (38) that these ultra-relativistic species do not contribute to the (F^) terms which source the dilaton fields' 
evolution. 

The only uncharged fundamental field is the neutrino, which we will assume to be so light that it remains relativistic 
at experimental and cosmological temperatures. Such neutrinos make only a very small contribution to the dilaton 
source terms. The corollary of this is that amongst the known fundamental matter species, we are justified in assuming 
that only those with non-zero charge contribute to the right-hand sides of equations (27)- (29). 

The see-saw mechanism for neutrino mass-generation results in three light neutrinos, masses m^ ^ (corresponding 

to the currently observed particles) and three very heavy neutrinos, masses M^' « a few tens of GeV. The light 
neutrinos are formed mostly from the weakly interacting left-handed components, whilst the heavy neutrinos are 
primarily composed of the weak singlet right-handed particles. Such heavy neutrinos would therefore interact with 
the weak bosons much more weakly that other massive particle species. To zeroth order in m^^ /M^^^ <C 1, we can 
assume these heavy neutrinos to be non-interacting with either electromagnetic or weakly interacting particles. Their 
contribution to the (F^) terms will be negligible compared to that of the other matter species. 



B. The Top Quark 



With a mass of about 180 GeV, the top quark is the only fermionic species present in the Standard Model that 
does not fall into the <C 1 category. Whilst in principle its contribution to the dilaton terms could be calculated, 
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it would be a difficult procedure which would require accounting for gauge boson self-interactions, Higgs boson 
fluctuations, and non-perturbative effects - as well as having to do QCD calculations. However, the results would 
actually have very little bearing upon most cosmological tests of this theory, since the background energy density of 
top quarks is entirely negligible when compared to that of all the other forms of matter. 



C. 'Light' Charged Matter 



We have just argued that, amongst the Standard Model particle species, the only non-negligible contributions to 
the dilaton source terms will come from particles with mass m and charge Q 0, which are light compared to Mw 
We must also require the particles to be relativistic, allowing us to restrict to the low-energy broken symmetric phase 
of the KM clcctroweak theories. This is appropriate for energies well below the lOOGeV level. At these energies, 
perturbation theory will be appropriate and non-abelian effects will be of sub-leading order. In this case then, the 
results of Section II should be valid. 

Writing the (W^y • W'') and {Y^i, ■ Y^^") quantities in terms of the low-energy physical fields we find: 



■ Fw±) + sin2 0w {-\F^m) + cos^ (-3^!) + sinOw cosOw {-^F^, 

0{d^ew) + 0{gw,gY) 



(42) 



(-iV-r'"')= cos'ew{-\F^m) + ^^'^^Ow{-\Fl)-sinewcosew{-\F,m-Fz) ■ i + o{d^ew) m 



The O (d^dw) terms produce only negligible corrections. The O {gw,9Y) symbol represents the additional terms 
that arise from gauge boson self interactions; these only contribute at sub-leading order. 

Thus for a non-relativistic particle species, tp^, of mass <C M^y, charge Qi 7^ 0, and weak isospin 
t^i we see that only (— l-Fem) ^'^^^ contribute to (42, 43) at leading order. The ( — (- 

^ ^ Fern ' Fz) tcrms are suppressed by relative factors of order of QTpr^^^Tg^j^ {t':M ^ 2Q sin^ 9w 

Q e^cosBw 1^ (^'^^ ^ ^'^') I'sspectively. The leading-order contributions of such a 'light' matter species to 

(42) and (43) therefore reduce to 



2" em 

and 



3 \Fw±\) and 



2^ Pi 



(44) 



where 5m is defined as the electromagnetic mass correction when the electric charge is some value e = const. 



D. Dark Matter 



There is a great deal of evidence from cosmological and astronomical data for the existence of dark matter which 
contributes about 27% of the gravitating mass density of the universe. Such matter must be non-baryonic, electro- 
magnetically non-interacting and non-relativistic ('cold'). It appears that the Standard Model lacks any suitable dark 
matter candidates. Weakly Interacting Massive Particles (WIMPS), such as the lightest putative supersymmetric 
partners in MSSM, are one possible class of candidates for dark matter. The masses of these particles tend to be 
of the order of a few tens of GeVs and so fall into the m « Mw category. Locally, if dynamically virialised, they 
will have keV energies which may allow them to be detected in underground nuclear recoil experiments. They are 
necessarily uncharged ( Q 7^ 0) but they can however interact weakly and as such contribute to the ( — j^J) and 
(~3 l-^W'itl) terms. Thus the leading-order contribution from these WIMPs to (42, 43) would be given by: 



/_iW •W'^N 



(-iK^y^n . =^tan2eiy^^ 

\ 4 f*" / wimp ^ 



(45) 



m 



wimp 



Ml 



^wimp ' 
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where we have defined Tw and Tz to be WIMP-model dependent 'structure' functions. These encode precisely how 
the WIMP's C parameters depend on the gauge boson masses. We expect \Ty^\ and \Tz\ to be <C 1. Not much more 
can be said about their structure as functions of the dilaton field without the aid of a microscopic model for dark 
matter. We shall therefore leave this discussion for a separate work. It should be noted also that, if Tw and Tz are 
no more than an order of magnitude or so smaller than the C,^.^ parameter of baryonic matter, then WIMP-based dark 
matter could well be the dominant factor in the cosmological evolution of the dilaton fields because the cosmological 
dark matter density in WIMPs is an order of magnitude greater than that of baryonic matter. 

V. DISCUSSION 

In this paper we have derived the form of the effective dilaton field equations for an electroweak theory with 
varying couplings in the presence of a background matter density. It is this form that is of most use for experimental, 
observational and cosmological tests of this theory. We took particular care in the identification of the physical gauge 
fields. In previous work on BSBM theory, the dilaton dependence of the soiircc term was incorrectly determined 
as a result of the auxiliary fields mistakenly taken to be the physically propagating modes. This error resulted in the 
statement: 

Where /^^ ~ 2(?[^a,^] was the auxiliary field strength, and the auxiliary photon field. 5m is defined as the 
electromagnetic mass correction when the electric charge is some value e = const. This lead to the conclusion that 
the leading-order dilaton field dependence went like: ^ = e~^'^. However, it is clear from equations (39)-(40) that in 
reality the leading-order dilaton dependence should be Ij = e+^*' and equation (46) should read: 

(-'-'u.r) = '^'-4p = e^^CP (47) 
\ e"' 4 / e'^ m'' 

Indeed equation (46) is problematic because, in the limit of zero electric charge, matter decouples from the photon, 

and so cannot possibility contribute to the F,^,^ term. This situation corresponds to ip ^ — oo. IIowc;ver, in this limit, 
the term as given by the right-hand side of equation (46) grows infinity large. Equation (47) shows the correct 
behaviour, and vanishes, as it should, it this limit. 

We have also seen that if the gauge bosons become massive then, for particles which are much lighter than the 
gauge boson in question, the leading-order dilaton dependence of the F^ term changes from g"^ {g is the physical 

gauge coupling) to — ~ ■\. For particles that are much heavier than the gauge boson in question, provided 

gauge 9 

perturbation theory is still valid at energies of the order of the particle's mass, the leading-order dilaton dependence 
remains as g-^. We also briefly discussed the complication of mparticle ~ Mgauge, whereby the leading-order dilaton 
dependence will be highly susceptible to the details of the matter model in question, and noted that this effect might 
be important in cosmology if dark matter is weakly interacting. 

In all of this analysis we assumed both that perturbation theory holds and that any non-abelian effects are negligible. 
Whilst this is true for electroweak theory at energies well below the Higgs boson mass, it will not be true for QCD. If 
we were to construct a BSBM-likc varying asirong = fisirong i'^^ equivalently, varying Aqcd) theory, we would expect 
the leading-order dilaton dependence of the -Fstrong term to come from a complicated function of g strong- Evaluating 
this function would, at the very least, require us to be able to predict quark and nucleon masses accurately via a 
QCD calculation which is not yet possible. In the absence of such calculations, varying Aqcd theories will be difficult 
to test accurately or make use of in the early universe except at the very highest energies. Fortunately, we do not 
encounter these problems in the electroweak KM theories. 

We conclude with a statement of the effective dilaton field equations we have derived. In the KM-I theory these 
are: 

where the sum is over all charged matter species (with m <C Mw)- The first term is identical to the properly evaluated 
source term in BSBM. Hence the effective KM-I theory differs (at low energies and densities) from BSBM only in the 

putative WIMP matter contribution. If we transform ip — > — then we can read off the explicit correspondence with 
the the studies in refs. [28, 29, 31, 34, 35, 36, 37, 38]. When Q/uj < we obtain slow logarithmic growth of the fine 



^1 



w 



Mw {vY 



■ sec 6w^z 



wxmp 

Mz {'pf 



(48) 
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structure 'constant' during the dust dominated era, as oc ln[ln(t + to)], constant but constant-(/? behaviour during 
the radiation and dark-energy dominated eras. \i Q/uj > then the solutions predict a much stronger evolution of a 
with time that is difficult to reconcile with the observational constraints. The sign of Q/lo is controlled by the sign of 
C, e [—1, 1] and of w. Positive C, corresponding to 'normal' matter dominated by the electrostatic contributions, and 
negative C, corresponding to matter (like superconducting cosmic strings, which have C, = —1) that is dominated by 
magnetic energy. Positive a; corresponds to a positive kinetic contribution to the energy by the y> field while negative 
(jj indicates that it is a ghost field, as in ref. [45]. 

In the KM-II theory the effective dilaton equations are: 



\ 



2 2^ «er» (y, X) 

— sm fviA — 



2 



^1 



ml 



w 



Mw {vY 



Mz {V, X) 



Pu 



— cos tiw [tp, X) = 

W2 Q!em 



CiPi H e^^ tan^ 9w^z 



m 



wimp 



U)2 



Mz{tp,x) 



(49) 



(50) 



A similar menu of possibilities exists for the sign of the leading term on the right-hand side of the Dip and equations 

as was the case for the KM-I theory discussed above. In the absence of the permitted dark matter contributions, this 
two-dilaton theory will only reduce to BSBM when uj2 sin^ 9w = cos^ 0w In all other cases 9w will vary and lead 
to an evolution of aem that is different from that of BSBM theory. Further cosmological consequences of these results 
will be explored elsewhere. 
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